The gap symmetry is marginal in high-T c superconductivity: Isotropic s-wave pairing 
leads to either s or d-wave gap depending on the Coulomb pseudopotential 
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We point out that in two-dimensional systems like high-T c superconductors, an isotropic s-wave 
electron-phonon interaction in the momentum decoupling regime leads to either s or d-wave super- 
conductivity depending on the conventional values of the Coulomb pseudopotential. It is therefore 
possible that Y Ba^CuzOi is an anisotropic d-wave and Bi^Sr-zCaCu-zOs an anisotropic s-wave su- 
perconductor both having the same pairing interaction and similar T c 's. A transition from s-wave 
to d-wave superconductivity may also be possible by adjusting the doping. 
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The symmetry of the order parameter is in the center 
of the debate on the origin of high-T c superconductiv- 
ity 0]. It is generally believed that the symmetry of 
the gap should also indicate the symmetry of the pair- 
ing mechanism. Studying this symmetry should allow 
one to determine whether cuprates are conventional s- 
wave superconductors with an interaction that could be 
mediated by phonons, or whether they are unconven- 
tional d-wave superconductors in which case the pair- 
ing interaction should be due to spin fluctuations |^] 
or to some other exotic mechanism. In this spirit one 
expects that BiiSriCaCu20% and Y BaiCu-sOi should 
both have the same symmetry, and their high critical 
temperatures should be intimately related to the sym- 
metry of the gap. 

Qualitative aspects of anisotropy in BiiSriCaCuiO^, 
like the correlation between the gap, the electronic den- 
sity of states (DOS) and the visibility of the dip above 
the gap anisotropies Q| , the strong enhancement of the 
anisotropy with temperature p|, the asymmetry of the 
SIN tunnel spectra where the dip is seen only at negative 
sample bias 11 and also the possible presence of a sec- 
ondary gap maximum or "hump" in the T — X direction 
0, have been shown to indicate that the pairing inter- 
action is isotropic s-wave in the momentum decoupling 
regime M . Momentum decoupling (MD) arises when the 
characteristic momenta exchanged during the pairing in- 
teraction are small compared to the characteristic mo- 
menta of the electronic density of states (DOS) varia- 
tions over the Brillouin zone. In that case the Eliash- 
berg equations in the momentum space tend to decou- 
ple and have different couplings in different regions of 
the Fermi surface. These couplings are proportional to 
the local DOS and therefore anisotropies of supercon- 
ductivity reflect the DOS anisotropies. Our analysis 
of the B%2Sr2CaCu20s phenomenology in this context 
points to an anisotropic s-wave gap [Q. However, vari- 
ous phase sensitive S and node sensitive ^| experiments 
on YBC12CU3O7 report evidence of a sign reversal of the 
order parameter supporting d- waves ]l0| , although ex- 
perimental contradictions persist p^] , |12| | . 

We will show in this letter that, contrary to the general 



belief, the symmetry of the order parameter does not rep- 
resents a crucial criterion for the symmetry of the pairing 
interaction, and is an energetically marginal parameter. 
We will see in particular that an isotropic s-wave inter- 
action (that could be mediated by phonons) in the mo- 
mentum decoupling regime can lead to either s-wave or 
d-wave superconductivity depending on parameters that 
are marginal for the pairing mechanism. These parame- 
ters are the magnitude of the Coulomb pseudopotential 
\x* and the relative importance of the characteristic mo- 
menta of the variations of the Coulomb pseudopotential 
fi* compared to the characteristic momenta exchanged 
during the pairing interaction. It is for example possi- 
ble that Bi2Sr2CaCu20s is an anisotropic s-wave and 
Y Ba2Cu-sOi an anisotropic d-wave superconductor both 
having the same attractive isotropic s-wave interaction in 
the momentum decoupling regime. In addition, it may be 
possible to switch from s-wave to d-wave superconductiv- 
ity and vice versa by doping since \i* depends sensitively 
on the doping. Of course it will also be clear that a 
d-wave gap symmetry does not imply a spin fluctuation 
pairing mechanism Q . 

Energetic irrelevance of the symmetry of the order pa- 
rameter has also been reported in a model proposed to de- 
scribe the "spin gap" in underdoped cuprates This 
property is also seen in another model in which the Fermi 
surface is divided in three independent pieces and the in- 
teraction has a low energy cut-off E3]. The possibility 
that phonons may help d-wave superconductivity which 
is primarily due to spin fluctuations has been reported in 
Ref. [|l5| , while more recently Nazarenko and Dagotto |l6} | 
proposed a specific Holstein model with nearest neighbor 
attraction which leads to d-wave superconductivity the 
spin fluctuations playing a secondary role. 

The Coulomb pseudopotential /1* represents the effec- 
tive repulsion between the paired electrons and is 
present in all conventional superconductors, acting essen- 
tially as an effective negative coupling (l9) and reducing 
slightly the isotope effect exponent. In a two dimensional 
system in the MD regime fi* becomes a crucial parame- 
ter having unexpected implications on the gap symmetry. 
We performed simulations on a BCS model, in which case 
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the gap is given by (we use the notations of 0) 

Where V'(g) = V A (q) + n*(q) is the sum of the attractive 
interaction, due for example to electron-phonon coupling 
V A , and the effective Coulomb repulsion /x*. We solve 
this equation iteratively on a grid 120 x 120 in the zero 
temperature regime. 

We consider an isotropic s-wave electron-phonon cou- 
pling having at small momenta a Lorentzian behavior 
as a function of the norm of the exchanged momen- 
tum V A (q) = -9 2 (1 + \q\ 2 /Q 2 c (e - ph))" 1 . In this 
spectrum the electron-phonon scattering is dominated 
by the processes which transfer a momentum smaller 
than Q c (e — ph). In our analysis Q c (e — ph) is the 
relevant parameter and the particular shape of the in- 
teraction is irrelevant. As for the repulsive interaction, 
we first consider the simple case of a momentum inde- 
pendent Coulomb pseudopotential (x*(q) = Hq (hard- 
core-like yet finite repulsion). Band structure effects 
are completely marginal for our discussion here, and 
for clarity we show results corresponding to the simple 
nearest neighbor tight binding dispersion at half-filling 
£g = — 2t[cos(k x ) + cos(fcy)] (the lattice spacing is taken 
equal to unity). All our discussion is unaltered by the 
inclusion of next nearest neighbor or other terms in the 
dispersion relation Q . 

We show in figure 1 some of the calculated momentum 
dependent gap functions on the Fermi surface, for differ- 
ent values of /Xq in both the d-wave and s-wave channel. 
When k varies from to n we cover a quadrant of the 
Fermi surface Here we take Q c (e — ph) » 7r/12 which 
places us in the momentum decoupling regime and is rel- 
evant for BiiSriCaCuiO% [Q. The momentum depen- 
dence of the gap in the s-wave channel at zero repulsion 
(upper full line) is precisely due to momentum decou- 
pling. 

When we introduce repulsion in the s-wave channel, 
the gap is reduced by a constant amount in all directions 
(triple-dot-dashed line) resulting therefore in an effective 
enhancement of the anisotropy. At a critical value of the 
repulsion which in the case considered in Fig. 1 is on the 
order /x*/g 2 ~ 0.21 (dot-dashed line), the gap becomes 
almost zero in the (71", w) direction, and we have a dis- 
continuous transition to a new gap symmetry structure 
with two nodes in the quadrant shown in Fig. 1 (dotted 
and dashed lines). In this new state the gap becomes 
independent on the magnitude of fx* and the areas of the 
Fermi surface at which the gap is positive are equal to 
the areas in which the gap is negative. We will see in 
the following that this new state is not the physical state 
occupied by the system because at these fx* the d-wave 
solution is energetically more favorable. 



The d-wave solution is also characterized by indepen- 
dence of the gap on the magnitude of fx*. The origin of 
the repulsion independence of the gap lies on the momen- 
tum independence of p* , on the equality of the Fermi sur- 
face areas with positive and negative gap, and on the two 
dimensional character of the system considered here. We 
can see this analytically by considering a circular Fermi 
surface. If the gap on the Fermi surface has the form 
= A cos(n</>) where n is the number of nodes and 
cj> is the polar angle in the usual polar coordinates, then 
the momentum independent repulsive contribution to the 
gap function becomes proportional to the integral |l7j 

Jo J A + A 2 cos 2 (n</>) 

which is identically zero. Therefore a gap changing sign 
periodically on a two dimensional system eliminates the 
effect of any momentum independent repulsion. 

We obtain in this way a qualitative understanding of 
the discontinuous nature of the transition in the s-wave 
channel in Fig. 1 (from dot-dashed live at fx* / g 2 = 0.21 
to dashed line at li* /g 2 = 0.22). We also remark on the 
following points that might generally apply to any two- 
dimensional system: First, it is impossible to get a zero 
gap somewhere on the Fermi surface without a variation 
of sign of the order parameter. Secondly, the system can 
only have nodes with a finite slope and therefore if there 
are nodes a linear T-dependence of the penetration depth 
at low temperatures is plausible Finally, when the 
gap changes sign, the area of the Fermi surface in which 
it is positive is equal to the area in which it is negative. 

In realistic situations fx* is momentum dependent (the 
repulsion is not hard-core-like) To introduce a smooth 
momentum cut-off for the variations of fx* we consider 
a structure analogous to that of the pairing interaction 
written at small q as fi *(q) = Li* Q (l + \q\ 2 /Q 2 c (Cb)) , in 
which case Q c (Cb) represents the characteristic range of 
the exchanged momenta (smooth cut-off) in the repulsive 
interaction. We considered various different momentum 
structures of p* and we checked that the following discus- 
sion remains valid with only slight quantitative modifi- 
cations pc| . The important parameter is the momentum 
cut-off in the repulsive interaction Q c (Cb) compared to 
that in the attractive pairing interaction Q c (e — ph). We 
show in Fig. 2 the s-wave and d-wave gap solutions for 
Q c (Cb) — 7r/4 and Q c (e — ph) = n/12. Now the d-wave 
solution is not repulsion independent but it still appears 
less sensitive to the repulsion than the s-wave solution. 
It is clear from Fig. 2 (and also Fig. 1) that at some crit- 
ical repulsion the absolute d-wave gap will become larger 
in average than the s-wave gap, and this will become the 
energetically favorable state. 

To find out which of the s-wave and d-wave solutions is 
the physical state of the system, one has to calculate the 
free energy gain due to the superconducting transition 
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pl| , the physical solution being that with the higher ab- 
solute free energy gain (the lower free energy) . Of course 
the solution with the higher free energy gain is that for 
which the integral of the absolute value of the gap on 
the Fermi surface is higher and this is also the solution 
with the higher T c . In figure 3 we show the evolution 
of the free energy gain (absolute value of the free energy 
due to the superconducting transition in arbitrary units) 
as a function of p^jg 2 in the case Q c (Cb) = 7r/4 and 
Q c (e—ph) = it/12. The s-wave solution is favorable when 
/1q/<7 2 is small but when this ratio takes values larger than 
a critical value of the order of ^o/g 2 ~ 0.15 then the d- 
wave solution becomes more favorable. If we were able to 
control fi* for example by doping then we could induce 
transitions from s-wave to d-wave superconductivity and 
vice versa. 

Notice that such s-d transitions occur in our system 
at conventional values of the Coulomb pseudopotential 
already seen in low-T c superconductors |l8| . The critical 

is very dependent on the ratio Q c (e—ph)/Q c (Cb), We 
studied the evolution of the free energy gain in the s and 
d-wave channels as a function of fig/g 2 for different val- 
ues of Q c {e — ph)/Q c (Cb). This allowed us to construct 
the phase diagram shown in Figure 4. We also made 
an analogous study in the case of different momentum 
structures of /i* having an effective momentum cut-off 
and the resulting phase diagrams are similar to that of 
Fig. 4 When for a given Q c (e - ph)/Q c (Cb) the 

ratio /J,* j g 2 is larger than the critical value showed in 
Fig. 4, the gap is d-wave. It is clear that the smaller 
the ratio Q c (e —ph)/Q c (Cb) the smaller is the value of 
the critical repulsion while for Q c (e —ph) of the same or- 
der as Q c (Cb) the d-wave solution is impossible and this 
is probably the case in low-T c metallic superconductors. 
Since in the MD regime Q c (e — ph) <C kp, we can also 
have Q c {e — ph) <C Q c {Cb), in which case s and d-wave 
superconductivity are energetically close and both states 
are physically acceptable depending on the precise value 
and structure of /Li* . We remark that if in our scheme we 
want to have both d-wave and high-T c , then it is better 
to be well in the MD regime and have a small critical /i* , 
since large /i* reduces T c . Notice that the evolution of T c 
follows essentially that of the free energy gain and from 
figures 2 and 3 one can conclude that the negative effect 
of /i* on T c is smaller in the case of d-waves than in the 
case of s- waves [ po| . 

It is important to notice that the d-wave state we ob- 
tain is anisotropic and its anisotropies are driven by the 
DOS anisotropies in the same way as for the s-wave state 
0. This is clear in Figures 1 and 2 where the d-wave 
solution away from the (7r/2,7r/2) direction has exactly 
the form of the s-wave solution. When we are in the MD 
regime we always have different couplings in different re- 
gions of the Fermi surface whatever the symmetry of the 
gap and a large part of our discussion of ARPES and 
tunneling in B^'S^C'dC^Os M remains valid even for 



the d-wave state. In particular all our discussion on the 
origin of the dip structure |^2| , the correlation of gap, dip 
and DOS anisotropies, and the asymmetry of SIN tunnel 
spectra (?]] is also true in the MD d-wave. We remark 
that in the case of a d-wave state obtained by spin fluc- 
tuations it was impossible to reproduce the dip above the 
gap by strong coupling calculations p| . 

We are grateful to Prof. L. Pietronero for stimulating 
discussions and encouragement. 
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Figure 1: The gap as a function of k x (or k y ) for a 
momentum independent repulsion fi* and Q c (e — ph) = 
7r/12. We consider the s-wave channel and fi* = 
(upper full line), /i* y ' g 2 — 0.05 (triple-dot-dashcd line), 
li* /g 2 = 0.21 (dot-dashed line), {i* j g 2 = 0.22 (dotted 
line), n* I ' g 2 — 0.30 (dashed line) and the d-wave channel 
for the same n*/g 2 ratios (lower full line). The d-wave 
solution is independent of /z* . 

Figure 2: The s-wave gap (full lines) and the d-wave 
gap (dashed lines) as a function of k x (or k y ) on a quad- 
rant of the Fermi surface, for three different values of 
the repulsion ' g 2 = 0.05, 0.15 and 0.25, in the case 
Q c (e - ph) = 7r/12 and Q c (Cb) = n/4. The smaller the 
ratio Ho/g 2 the larger is the absolute value of the gap in 
both the s and d-wave solutions. 

Figure 3: The absolute value of the free energy gain 
due to the superconducting transition as a function of 
the ratio p%jg 2 for Q c (e - ph) = it/ 12 and Q c (Cb) = 
7r/4. The dotted line (triangles) corresponds to the d- 
wave solution and the full line (circles) to the s-wave 
solution. The physical solution is that with the higher 
absolute free energy gain (lower free energy). 

Figure 4: Phase diagram. The critical repulsion 
/ijj I g 2 for the transition from s-wave to d-wave supercon- 
ductivity as a function of the ratio Q c (e — ph)/Q c (Cb). 
The upper region of the graph corresponds to the d-wave 
state. 
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